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ABSTRACT 

We examine the clustering properties of low-power radio galaxies at redshift 0.4 < 
z < 0.8, using data from the 2SLAQ Luminous Red Galaxy (LRG) survey, and find 
that radio-detected LRGs (with typical optical luminosities of 3 — 5L* and 1.4 GHz 
radio powers in the range 10^'' to 10^^ WHz^^) are significantly more clustered than a 
matched population of radio-quiet (< 10^^ WHz^^) LRGs with the same distribution 
in optical luminosity and colour. 

The measured scale length of the two-point cross-correlation function between 
the full LRG sample and the radio-detected LRGs is 9.57±0.50/i~^ Mpc, com- 
pared to 8.47±0.27ft-^^ Mpc for the matched sample of radio-quiet LRGs; while 
the implied scale length of the auto-correlation function, tq, is 12.3±1.2 Mpc 
and 9.02±0.52/i~^ Mpc for the radio-detected and radio-quiet samples respectively. 
We further interpret our clustering measurements in the halo model framework and 
demonstrate that the radio-detected LRGs have typical halo masses of 10.1±1.4 
xlO^^h-^MQ and bias of 2.96±0.17, compared to 6.44±0.32 xW^^h'^MQ and 
2.49±0.02 for the radio-quiet sample. A model in which the radio-detected LRGs 
are almost all central galaxies within haloes provides the best fit to the measured 
clustering, and we estimate that at least 30% of all 2SLAQ LRGs with the same 
clustering amplitude as the radio-detected LRGs are currently radio-loud. 

Our results imply that radio-detected galaxies in the 2SLAQ LRG sample typ- 
ically occupy more massive haloes than other LRGS of the same optical luminosity, 
so the probability of finding a radio-loud AGN in a massive galaxy at z ~ 0.55 is 
influenced by the halo mass and/or cluster environment in addition to the well-known 
dependence on optical luminosity. If we model the radio-detected fraction of LRGs, 
Frad, as a function of halo mass M, then the data are well-fitted by a power law of the 
form Fiad c)c m"-^^*'^-^'^. The observed relationship between radio emission and cluster- 
ing strength could plausibly arise either through a higher fuelling rate of gas onto the 
central black holes of galaxies in the most massive haloes (producing more powerful 
radio jets) or through the presence of a denser IGM (which would provide a more 
efficient working surface for the jets, thus boosting their observed radio luminosity). 
Further work is needed to determine which of these effects is dominant. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: active 
- radio continuum: galaxies - galaxies: evolution - galaxies: elliptical and lenticular, 
cD - cosmology: large-scale structure of Universe 



1 INTRODUCTION 

1.1 Radio galaxies and their environment 

It has long been known that the hosts of powerful radio-loud 
AGN are massive early-type galaxies, and that the probabil- 
ity of such a galaxy hosting a radio source increases rapidly 
with optical luminosity and stellar mass (Auriemma et al. 
1977; Best et al. 2005a; Mauch & Sadler 2007). What re- 



mains less clear is the role (if any) of a galaxy's environment 
in determining whether it hosts a radio source. 

Studies carried out in the 1980s hinted at significant 
differences between the clustering properties of the powerful 
radio galaxies and the less powerful FRI sources. 



^ Fanaroff & Riley (1974) divided radio galaxies into two classes 
based on their observed radio morphology. They found a correla- 
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Based on an imaging study of 43 radio galaxies at 2 < 
0.3, Heckman et al. (1986) found that the most powerful 
radio sources (with radio luminosities above lO^'' "' W Hz~^ 
at 408 MHz) lay in regions of below-average galaxy den- 
sity and often showed a disturbed morphology suggestive 
of a recent interaction with a gas-rich companion. In con- 
trast, less powerful radio sources appeared to be associated 
with morphologically-normal eaxly-type galaxies in regions 
of high local galaxy density. Similar results were obtained 
by Prestage & Peacock (1988), who used an angular cross 
correlation technique to study the clustering environment 
of a sample of about 200 nearby {z < 0.15) radio galax- 
ies. They found that FRI radio galaxies lay in regions of 
significantly enhanced galaxy density, whereas the cluster- 
ing environment of FR II sources was similar to that of the 
overall population of 'normal' elliptical galaxies. 

The picture changed significantly in the mid 1990s with 
the work of Lodlow & Owou (1996), who discovered that 
the division in radio power between FRI and FRII radio 
sources was a strong function of the optical luminosity of 
the host galaxy. As a result, FRII radio sources are gen- 
erally hosted by less optically luminous (and less massive) 
galaxies than FRI sources of similar radio power. Since more 
massive galaxies also tend to be more strongly clustered, 
this effect needs to be taken carefully into account when 
analysing the clustering properties of powerful radio sources. 
Ledlow & Owen (1996) measured the bivariate radio lumi- 
nosity function (RLF) of early type galaxies in rich clusters, 
and found no statistically significant difference between the 
RLFs of galaxies in rich clusters and in the field. Their re- 
sults suggested that the local environment plays little or no 
role in determining whether an early-type galaxy hosts a 
radio-loud AGN, and that the optical luminosity and other 
properties of the host galaxy are by far the most important 
parameters affecting radio source formation and evolution. 

1.2 Measurements of radio— source clustering 

The advent of large galaxy rodsliift surveys like the 2dF- 
GRS and SDSS (Colless et al. 2001; York et al. 2000), com- 
bined with "all-sky" radio continuum surveys like NVSS 
and SUMSS (Condon et al. 1998; Bock et al. 1999) made it 
possible to assemble samples of thousands of objects with 
which to carry out statistical analyses of radio galaxies in 
the local universe (Best et al. 2005b; Mauch & Sadler 2007). 
The NVSS and SUMSS source catalogues are large and uni- 
form enough that the characteristic imprint of large-scale 
structure can easily be seen in the angular correlation func- 
tion (Blake & Wall 2002; Blake et al. 2004). Convolving the 
angular clustering amplitude in these surveys with a charac- 
teristic redshift distribution N{z) suggests that the present- 
day clustering length ro of radio galaxies is in the range of 
7-10/i^^Mpc, corresponding to a clustering strength sim- 
ilar to optically-luminous elliptical galaxies in moderately 
rich environments (see e.g. Overzier et al. 2003). 

Recently, Best et al. (2007) have revisited the question 
of radio-source clustering using data sets much larger than 

tion between morphology and radio luminosity, with less luminous 
(FRI) sources having a jet- like appearance and more luminous 
(FRII) sources having edge- brightened radio hotspots. 



those available to Ledlow & Owen (1996). Using a sample 
of 625 nearby galaxy groups and clusters selected from the 
SDSS, they show that the brightest galaxies in groups and 
clusters (BCGs) are more likely to host a radio-loud AGN 
than other galaxies of the same stellar mass. The probability 
is increased by up to a factor of two for the most massive 
galaxies (with stellar mass ~ 5 x 10^^ Mq), and by over an 
order of magnitude for galaxies with stellar masses below 
10^^ M0. This enhanced likehhood of radio-loud AGN ac- 
tivity was only seen in the innermost regions of a group or 
cluster (i.e. within 0.2r2oo, where rgoo is the Virial radius 
of the cluster). Best et al. (2007) argue that the radio prop- 
erties of both BCGs and non-BCGs can be explained if the 
radio emission is mainly fuelled by cooling gas from an X 
ray halo surrounding the galaxy. It therefore appears that 
although the radio properties of most galaxies in the local 
universe are largely unaffected by their environment, this is 
not true for massive galaxies located in the central regions 
of clusters. 



1.3 The 2SLAQ LRG radio sample at z ~ 0.55 

At higher redshift, Sadler et al. (2007) recently combined 

data from the 2SLAQ Luminous Red Galaxy (LRG) redshift 
survey (Cannon et al. 2006) and the NVSS and FIRST radio 
surveys (Condon et al. 1998; Becker et al. 1995) to identify 
a volume-limited sample of 391 radio galaxies at redshift 
0.4 < z < 0.7. They measured the redshift-space correla- 
tion between the radio- detected 2SLAQ LRGs and the full 
LRG sample, and found that the 2SLAQ radio galaxies were 
more strongly clustered than the overall 2SLAQ LRG pop- 
ulation. Since the 2SLAQ radio galaxies as a class were also 
more optically luminous than the overall LRG sample it was 
unclear whether the increased clustering was a luminosity ef- 
fect, or represented a genuine difference in the environments 
of radio-loud and radio-quiet LRGs at « ~ 0.55. In the 
current paper, our goal is to answer this question by investi- 
gating the clustering properties of radio-loud 2SLAQ LRGs 
in more detail. 



1.4 Radio— galaxy duty cycles 

The inferred lifetimes of the radio sources associated with 
massive galaxies (typically 10^-10* yr; Parma et al. 1999) are 
significantly shorter than the ages of their parent galaxies, so 
it is generally assumed that all massive galaxies must cycle 
between radio-loud and radio-quiet phases over time. Feed- 
back mechanisms in which the hot intergalactic gas episod- 
ically cools to fuel a central AGN, and is then reheated by 
the ensuing radio jets (e.g. Binney & Tabor 2005, Ciotti & 
Ostriker 2007) provide a natural explanation for this pro- 
cess. 

For the 2SLAQ LRG sample, .Johnston et al. (2008) 
have shown that the stellar populations of radio-detected 
and radio-quiet galaxies are generally indistinguishable. 
This is consistent with a picture in which 'radio-mode' AGN 
feedback (Bower et al. 2006; Croton et al. 2006) regulates 
the star-formation rate in these massive galaxies, and all of 
them undergo radio-loud episodes when their central black 
hole is active and can power radio jets. If this is the case, 
then the fraction of galaxies which are detected as radio 
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Table 1. Summary of the three 2SLAQ LRG samples used in 
this clustering study. 



Set 



Properties 



N 



1 All 'Sample 8' LRGs with good-quality z 7009 

2 Radio-detected LRGs from Set 1 250 

3 Luminosity— matched set of LRGs from Set 1 2750 



sources above some radio power Pum simply represents the 
fraction of the radio-galaxy duty cycle for which a typical 
galaxy is a radio source at or above this level. 

Observationally this is complicated by the fact that the 
radio luminosity function of AGN is very broad, spanning 
at least six orders of magnitude (Mauch & Sadler 2007), 
and there is also a strong correlation between radio power 
and optical luminosity (Auriemma et al. 1977, Best et al. 
2005). The observed radio detection rate will therefore de- 
pend strongly on both radio power and galaxy luminosity. 
In this paper we consider only the very luminous early- 
type galaxies which comprise the 2SLAQ LRG sample, so 
that the range in optical luminosity is small. We also use 
the term 'radio-detected' to refer to galaxies whose 1.4 GHz 
flux density is higher than the 1-2 mjy detection limit of the 
FIRST and NVSS radio surveys. Since the 2SLAQ LRG ra- 
dio sample is close to volume-limited (see Figure 7 of Sadler 
et al. 2007) , this translates to a limiting radio luminosity of 
~ 102*-2WHz-\ 

Throughout this paper, we assume a flat A-dominated 
cosmology with f2m=0.27, Ho=70 km s~^Mpc~^, and 
(78=0.8 unless otherwise stated. 



2 THE 2SLAQ LRG DATA SAMPLES 

In this paper we will consider the clustering of three data 
sets, as summarised in Table [T] All three samples are drawn 
from the full 2SLAQ LRG spectroscopic survey (Cannon 
et al. 2006). The first data set consists of all main sample 
('Sample 8') LRGs with high quality redshifts and 0.45 < 
z < 0.8 (see Cannon et al. 2006 for details). The second 
set is a sub-sample of the first, and includes only the LRGs 
which have been detected as radio sources (see Sadler et al. 
2007). Since we want to test whether radio-detected LRGs 
are more strongly clustered, we select a further sub-sample 
which has the same optical properties as the radio-detected 
sub-sample but does not contain any LRGs with detected 
radio emission. 

Figure [1] shows the distribution of the intrinsic lumi- 
nosity, redshift and colour of the whole LRG sample and 
the radio-detected subsample, where we have used the K-l-e 
corrections as detailed in Wake et al. (2006) to generate the 
intrinsic luminosity and colours. This figure shows that the 
radio sample is intrinsically more luminous, has a slightly 
higher typical redshift and is slightly redder. As a result of 
this, to generate a matched sample of LRGs without radio 
emission we must match for both optical luminosity and red- 
shift. To be complete we also match for colour, even though 
the difference in the distributions is marginal. We note that 
excluding this colour matching makes no difi'erence to our re- 
sults. We generate our sample by selecting the eleven LRGs 
from the whole sample that lie closest to each radio-detected 



LRG in colour-magnitude-redshift space. The final matched 
sample contains 2750 LRGs, which is the maximum that 
can be generated if we want to match the radio distribu- 
tions without a large number of repeats (~15%). 



3 THE TWO-POINT CROSS-CORRELATION 
FUNCTION 

Since the space density of the radio LRG sample is very 
low, measurements of its auto-correlation function would be 
dominated by shot noise, particularly on small scales. In 
order to reduce this noise we cross-correlate the radio and 
matched samples with the full LRGs sample, which has a 
much higher space density. 

The 2pt-cross-correlation function between two sets of 
objects a and 6, ^ab{r), is defined as a measurement of 
the excess probability above Poisson of finding an object 
a at a separation r from another object b. Here we wish to 
cross-correlate the radio and matched samples with the full 
LRG sample. We calculate this by comparing the number of 
pairs as a function of scale between the radio-detected (or 
matched) and full sample, with the number of pairs between 
the radio (or matched) and an unclustered (random) cata- 
logue, which covers the same volume as the full sample such 
that 

nRNRf(r) 



C(s) = 



(1) 



rif NRr{r) 

where Njif and Nrt are radio- full and radio-random pair 
counts respectively, and n/ and are the number of galax- 
ies in the full and random samples. 

When making this calculation for our samples we must 
take into account the effect of the completeness varying 
across our survey. We follow the procedure described in de- 
tail by Wake et al. (2008), by up- weighting LRGs in areas of 
low completeness and using a random catalogue that has a 
constant space density over the angular mask of the survey. 
We exclude from our calculations regions that have < 65% 
completeness or that are close to bright stars. 

We estimate the errors on our 2pt-cross-correlation 
function measurements using jack-knife re-sampling (Scran- 
ton et al. 2002; Zehavi et al. 2005). We split the 2SLAQ 
area into 74 equal area regions, minimising the noise on the 
covariance matrix whilst still removing regions larger than 
the scales we are interested in. We then repeatedly calculate 
each 2pt-function removing one area at a time to generate 
a full covariance matrix. Throughout we generate the pair 
counts using the KD-tree code in the NTROPY software 
package (Gardener et al. 2007). 

The peculiar velocities of galaxies generate an error 
in the distance measurement to a galaxy along the line of 
sight, which results in distortions to ^ known as redshift 
space distortions. To remove this efi'ect we can calculate the 
clustering perpendicular (r^) and parallel (tt) to the line-of- 
sight {^{rp,n)) and then integrate over the tt direction to 80 
/i~^Mpc to give the projected correlation function (w(rp)) 
such that 



Wp{rp) 



(2) 



This can be expressed in terms of the real space correlation 
function ^(r) (Davis & Peebles 1983) with 
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Figure 2. The redshift space 2pt-auto-correlation function for 
the full LRG sample (stars) and the redshift space 2pt-cross- 
correlation function for the radio sample (open circles) and a 
sample randomly selected from the whole sample to match the 
luminosity, colour and redshift distributions of the radio sample 
(filled circles). 

p oc 

Wp{rp) = 2 rdr^{r){r^ - rl)-^^\ (3) 

J rp 

If a power law of the form ^(r) — (ro/r) ' is assumed 
then equation [3] can be solved analytically (Davis & Pee- 
bles 1983). 



4 CLUSTERING PROPERTIES OF 

RADIO-DETECTED LRGS AT Z ~ 0.55 

Figures [2] and [3] shows the 2pt-auto-correlation function for 
the full sample, and the 2pt-cross-correlation function for the 
radio-detected and matched samples in redshift space and 



Figure 3. The projected 2pt-auto-correlation function for the 
full LRG sample (stars) and dark matter (triangles), and the 2pt- 
cross-correlation function for the radio sample (open circles) and 
a sample randomly selected from the whole sample to match the 
luminosity, colour and redshift distributions of the radio sample 
(filled circles). The lines show power— law fits on scales 0.5 < Vp < 
5Qh~^ Mpc with slope, 7, fixed at 1.83 corresponding to the best 
fit value for the full 2SLAQ LRG sample. 



projection respectively. In order to compare the clustering 
strengths of these samples we make fits to w{rp) assuming 
a power law form for ^(r) of (r/ro) ' using the analytical 
solution to equation [S] We use the full covariance matrices 
generated using the jack-knife resampling technique and fit 
over 0.5 < < 50h~^ Mpc. Values of ro and 7 for these 
fits are given in Table (2] Since the values of the slope, 7, 
are consistent between the three samples we refit using the 
7 from the fit to the full LRG sample. The values of ro from 
these fits are also given in Table[2]and the Ax^ distributions 
are shown in Figure U 
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Figure 4. Ax^ for the power law fits to the projected 2pt-cross- 
correlation functions for the radio sample (blue dashed) , a sample 
matching the luminosity, colour and redshift distributions of the 
radio sample (red solid), and the auto-correlation function of the 
full LRG sample (black dot-dashed). The slope, 7, is fixed at 
1.83 corresponding to the best fit value for the full sample. The 
horizontal dotted lines show the 68%, 90% and 99% confidence 
intervals. 



Table 2. Values of the power-law fits to the projected 2pt-auto- 
correlation function for the full sample and dark matter, and the 
2pt-cross-correlation function for the radio and matched samples 
in the range 0.5 < r < 50/i^^ Mpc. Errors are at the 68% confi- 
dence level. 

Sample ro (h"! Mpc) 7 x^ed ''0 (fixed 7) 

7.66l»;i« 1.83lg;°^ 7^661^ 

Matched 8.47tJ5'^^ l-78too5 8.3lt°'^^ 

Radio 9.57lgjJ l.7^+J^'^ 1.8 9.72+J^% 

Dark Matter ?,.2btl'°o2 l-81-o:o2 



Figure |4] and Table [2] show that the radio-detected 
LRGs are more clustered (with a significance of 97%) than 
other LRGs with similar optical luminosity. When consid- 
ering Figures 2, 3 and 4 and the fits given in Table [2j it 
is important to remember that the amplitude of the cross- 
correlation of the radio or luminosity matched samples will 
be lower than the auto-correlation for these samples, since 
they are cross-correlated with the full LRG sample which has 
a lower clustering strength. On large (> 1 Mpc) scales, where 
the clustering is determined by pairs of galaxies in separate 
dark matter haloes, one would expect the cross-correlation 
function to be the geometric mean of the auto-correlation 
of the two samples. On smaller scales where the clustering 
is dominated by pairs within haloes, the cross-correlation 
function will depend on the relative distribution of the two 
samples within haloes and would only be the mean if the 
distribution was identical. We can therefore use the cross- 
correlation on large scales to calculate the auto-correlation 
function, where 



Table 3. The scale length of 2pt-auto-correlation function (rp) 
for the three samples along with the inferred large scale bias and 
typical dark matter halo mass. Errors are at the 68% confidence 
level. 



Sample 


ro (h~^ Mpc) 


bias 


Mdh(1O"/j~'M0) 


Full 


7.66±0.17 


2.14±0.05 


3.1±0.2 


Matched 


9.02±0.52 


2.52±0.16 


5.7±1.1 


Radio 


12.3±1.2 


3.11±0.29 


10.6±1.9 



C2=C?2/Cl- (4) 

This gives values of ro for the auto-correlation function of 
12.3 ± 1.2 and 9.02 ± 0.52 h'^ Mpc for the radio-detected 
and matched samples respectively. 

We can estimate the large scale bias (6) for our three 
populations as & = {WyOai /wpDuY^^ where WpOM is gener- 
ated using the linear theory power spectrum as described 
in section [5] and the projected auto-correlation function 
for the matched and radio-detected samples is calculated 
using Equation 4. We define the large scale bias here as 
the weighted mean of [WpGai/wpOMY^'^ within 2 < r < 
50/1"^ Mpc. We can then use the dependence between the 
dark matter halo mass and bias estimated from the halo 
mass function (Sheth & Tormen 1999) to relate this bias to 
a typical halo mass. Table [3] gives the mean bias and halo 
mass for each of the three samples. 

Our finding that the radio-detected 2SLAQ luminous 
red galaxies (which have typical 1.4 GHz radio powers 
of 10^* to lO^^WHz-i) are significantly more strongly- 
clustered than LRGs of similar optical luminosity which 
are not detected as radio sources (and so are weaker than 
~ 10'^'*WHz~^) implies that the radio-galaxy duty cycle 
z ~ 0.55 is affected by at least one factor which is linked to 
the clustering environment, and we investigate this further 
in the next section. 



5 HALO MODELS 
5.1 Model parameters 

The halo model assumes that the galaxy clustering signal 
encodes information about the Halo Occupation Distribu- 
tion (HOD;how the galaxies populate dark matter haloes), 
in particular how the HOD depends on h alo mass (see e.g. 
Jing et al 1998, Ma & Fry 2000, Peacock & Smith 2000, 
Seljak 2000, Scoccimarro et al. 2001, Berlind & Weinberg 
2002). We have successfully applied this technique to the 
2pt-correlation function of the 2SLAQ LRGs (Wake et al. 
2008), and use the same techniques here to gain a further un- 
derstanding of how the radio galaxies are distributed within 
dark matter haloes. We give below a brief outline of our halo 
model description of the clustering, and refer the reader to 
Wake et al. (2008) for further details. 

In the halo model, every galaxy is associated with a halo 
and all haloes are 200 times the background density what- 
ever their mass M. Sufficiently massive haloes typically host 
more than one galaxy. The halo model we use distinguishes 
between the central galaxy in a halo and the others, which 
are usually called satellites. 
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Table 4. Details of the HOD fits to the auto-correlation function for the full sample and the cross-correlation function for the matched 
and radio samples. 



Selection 


Density 




Ml 


a 






Mcff 


Fsat 






(10i3?t-iMQ) 


(10i3?t-iMQ) 






(1013/i-iMq) 


(%) 


Full 


1.55 ± 0.02 


1.96 ± 0.02 


24.6 ± 1.1 


1.93 ± 0.13 


1.2 


2.15 ± 0.01 


4.56 ± 0.09 


5.2 ± 0.7 


Matched 


0.45 ± 0.01 


4.70 ± 0.09 


> 59.5 


> 0.85 


1.0 


2.49 ± 0.02 


6.44 ± 0.32 


0.9±1.9 


Radio 


0.15 ± 0.05 


9.65 ± 2.2 


> 80.0 


> 0.6 


0.82 


2.96 ± 0.17 


10.1 ± 1.4 


0.0 ± 4 
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Figure 5. The projected 2pt-auto-correlation function for the full LRG sample (stars) and the projected 2pt-cross-correlation function for 
the radio sample (open circles) and a sample randomly selected from the whole sample to match the luminosity, colour and z distributions 
of the radio sample (filled circles). The solid lines show the best fitting HODs on scales 0.32 < rp < 5Qh~^ Mpc. The dashed lines show 
the fit to the full sample, and the dotted line the fit achieved applying a power-law radio fraction to the matched HOD fit (see text for 
details). 



The fraction of haloes of mass M which host centrals is 
modelled as 



{Nc\M) = exp(-Af™i„/M). 



(5) 



Only haloes which host a central may host satellites. In such 
haloes, the number of satellites is drawn from a Poisson 
distribution with mean 

(Ns\M) = {M/M^r. (6) 
Thus, the mean number of galaxies in haloes of mass M is 
(iV|A/> = {N,\M)[1 + {Ns\M)], (7) 
and the predicted number density of galaxies is 



dMn{M) {N\M), 



(8) 



where n{M) is the halo mass function, for which we use the 
parametrisation given by Sheth & Tormen (1999). 

We further assume that the satellite galaxies in a halo 
trace an NFW profile (Navarro et al. 1996) around the 
halo centre, and that the haloes are biased tracers of the 
dark matter distribution. The halo bias depends on halo 
mass in a way that can be estimated directly from the halo 
mass function (Sheth & Tormen 1999) . With these assump- 



tions, the halo model for ^(r) is completely specified (e.g. 
Cooray & Sheth 2002). Our model for the real space 2pt- 
auto-correlation function is given in detail in Wake et al. 
(2008), where the mean number density of central-satellite 
pairs within haloes of mass M is n{M) {Nc\M) {Ns\M), and 
the mean number density of distinct satellite-satellite pairs 
is n(M) {Nc\M) {Ns\Mf/2. In this work we use the 2pt- 
cross-correlation function, where we are cross-correlating a 
sub-sample of LRGs (either the radio-detected or matched 
samples) with the full sample of LRGs. We model the cross- 
correlation function in the halo model framework analo- 
gously to the auto-correlation function with the mean num- 
ber density of central-satellite pairs. 



nc.(M) = n{M) {iVci|M) [(iV.|M) + {iVsi|M)], 



(9) 



and the mean number density of distinct satellite-satellite 
pairs. 



nss{M) = n{M) {Nci\M) (7V.|M) {7V.i|M), 



(10) 



where the pairs are between the full sample and sub-sample, 
and the terms with a subscript 1 are for the sub-sample and 
those without are for the full sample. 

Since in our definition of the HOD a halo must contain a 
central if it is to contain a satellite, whether for the full sam- 
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pie or sub-sample, then only haloes that have a central in the 
sub-sample will contribute to the cross-correlation function. 
This is why only the fraction of haloes which contain a cen- 
tral of the sub-sample {Nci\M) appears in equations |9] and 
[TOl not the equivalent fraction for the full sample {Nc\M). 

Our model for the real-space 2pt-cross-correlation func- 
tion is then 



where 

l+£,csi2{r) = 



1 + S,ssi2{r) 
and 

6hl2('") = 

with 

P2hl2{k) 



(r) + I + iss{r) + i2hi2{r) 



dM 



dM 



ncsu{M) p{r\M) 

nging2 M 

nsswjM) \{r\M) 
nging2 



dk lC'P2hi2{k) sinkr 
k 27r^ kr 



bgi{k)bg2{k)Pu^{k), 



(11) 

(12) 
(13) 

(14) 
(15) 



where 

bg{k) = 



dM^b{M) 



{N^\M) + {Ns\M)u{k\M) 



In the expressions above, p{r\M) is the density profile of 
haloes of mass M, \{r\M) denotes the convolution of two 
such profiles, u{k\M) is the Fourier transform of p{r\M) / M , 
and PLin(fc) denotes the linear theory power spectrum. In 
practise, we approximate hg(k) by its value feg at fc = 
(eauation ll6l) . All these quantities, along with the mass func- 
tion n(m) and halo bias factor b{M), are to be evaluated at 
the redshift of interest. We then calculate w{rp) from ^(r) 
using equation ((3]). 

5.2 Application to the 2SLAQ LRG data 

This model for the cross-correlation-function makes the as- 
sumption that the two galaxy samples occupy the same 
haloes, and that the satellite galaxies of both samples follow 
the same profile within the haloes. Whilst this is not nec- 
essarily valid for two independent galaxy samples, it should 
hold here as one sample is always a subset of the other. 

This form of the halo model makes the explicit assump- 
tion of a volume-limited sample of galaxies. The 2SLAQ 
LRG sample as a whole is magnitude-limited rather than 
volume-limited, but (as may be seen from Table 4 of 
Sadler et al. 2007) , the radio-detected LRGs are significantly 
brighter than the LRG sample as a whole. For this reason, 
Sadler et al. (2007) argue that the sample of radio-detected 
2SLAQ LRGs is close to volume-limited, with no strong cor- 
relation between absolute magnitude (or radio luminosity) 
and redshift. 

Nevertheless, since the use of a magnitude-limited sam- 
ple may have the effect of biasing any derived quantities we 
may wish to measure from the HOD (such as the typical 
halo mass or satellite fraction) we carried out some further 
investigations. In the Appendix, we investigate the conse- 
quences of applying the halo model to magnitude-limited 
galaxy samples using the latest semi-analytic galaxy forma- 
tion models (Font et al. 2008) applied to the Millennium 




^Ihalo (h-'M^) 



Figure 6. The mean number of LRGs per halo as a function of 
halo mass (top) and the mean number of LRGs per halo nor- 
malised by the number of haloes as a function of mass (bottom) 
for the best fitting HODs of the three samples. The full sample is 
shown by the black lines, the radio sample is shown by the blue 
lines and the luminosity, colour and redshift matched sample is 
shown by the red line. The total, central and satellite contribu- 
tions are shown by the solid, dashed and dotted lines respectively. 
Whilst the spike at large halo masses in the radio HOD is the best 
fit, the HOD parameters describing the satellites are very poorly 
fit and a model with no satellites represents a better fit to the 
data. 



simulation (Springel et al. 2005). We show that the halo 
model successfully recovers the effective halo mass to within 
5% and so we are confident in applying this model to the 
2SLAQ LRG samples. 

We can now fit this model to the measured clustering 
and find the best fitting HOD parameters Mmin, Mi, and a. 
We first make a fit to the full LRG sample, fitting both 
the measured 2pt-auto-correlation function and the mea- 
sured space density of 1.55 x 10~*/i^Mpc~'^ as was done by 
Wake et al. (2008). We then use the resulting best fit HOD 
when fitting the cross-correlation function of the matched 
and radio samples. We note that including the density for 
the radio sample in this fitting process produces a poor fit 
to the clustering. This is not surprising, since only a small 
fraction (3-4%) of the 2SLAQ LRGs are detected as radio 
sources and so their overall space density is very low. Our 
HOD parametrisation assumes that all haloes above a cer- 
tain mass contain at least one galaxy. Therefore the lower the 
galaxy density, the higher the typical mass and the higher 
the clustering. This is a good parametrisation for the whole 
galaxy population limited by optical luminosity, due to the 
correlation between central galaxy stellar mass and the halo 
mass. However, if the duty cycle of radio emission is such 
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that only a fraction of LRGs that could be radio loud actu- 
ally show emission at a given time, then the measured space 
density of the radio-detected objects is a product of both 
the parent population space density and the fractional duty 
cycle. 

To account for this, we could modify our HOD to in- 
clude a term determining the fraction of LRGs which are 
currently radio-emitting. Initially we chose a simpler ap- 
proach by not including the density and fitting the cluster- 
ing alone. We can then compare the density predicted by 
the best-fitting HOD to the measured value and determine 
the fraction of LRGs with the same clustering as the radio 
sources which are currently radio-loud. By doing this, we 
are assuming that the relationship between halo mass and 
the probability of a galaxy having detectable radio emission 
takes the same form that we use to relate halo mass and lu- 
minosity, multiplied by some fraction which is independent 
of halo mass. 



5.3 Results 

Figure [5] shows the best fitting HODs for the three samples 
with the best-fitting parameters given in Table |4] We can 
use the HODs to calculate some other useful quantities such 
as the average linear bias (bg), the effective halo mass (Mefi) 
and the satellite fraction (Fsat), where 



M, 



and 



Fsat 



dMn(M)b{M) <N > /ug 



dMMn{M) < N > /ug 



J dMn{M) <N,\M > 
J dMn{M) < N\M > ' 



(16) 



(17) 



(18) 



The values of these quantities for each sample are given in 
Tabled 

As we would expect from the relative clustering 
strength, we find a significantly higher minimum halo mass 
(Afmin) for the radio sample than the matched sample, con- 
firming that the radio galaxies are typically found in more 
massive haloes. This is further confirmed by the calculated 
values of Moh and bg. 

The parameters Mi and a describing the satellite pop- 
ulation are quite poorly constrained by these fits, with a 
large range of parameter space with high Mi or a having 
a very similar x^- This is because these samples have very 
low numbers of satellites, with the best fitting HODs for the 
matched and radio samples being consistent with having no 
satellite galaxies. 

The best-fitting HODs are plotted in Figure (6] which 
shows the number of galaxies as a function of halo mass 
(top) and the number of galaxies weighted by the number of 
haloes as a function of halo mass (bottom). Each plot shows 
the total galajcy distributions, as well as those of the cen- 
tral and satellite galaxies independently. These plots clearly 
show that for both the radio-detected and matched sam- 
ples, satellite galaxies are only present in the most massive 
haloes. Since there are so few haloes with these high masses, 
their total contribution is very low. This also explains why 
Ml and a are given as la lower limits in Table |4] When 



Ml becomes significantly large there is essentially no con- 
tribution to the clustering from satellites since there are so 
few haloes at these masses, so the fit is equally as good for 
any large value of Mi or a. In fact, fitting a model with no 
satellite galaxies produces a only marginally higher for 
both the radio-detected and matched samples, and so gives 
significantly better fits overall because of the larger number 
of degrees of freedom in the model with satellites. 



5.4 HOD models and the radio— galaxy duty cycle 

We derive a space density 0.15 ± 0.05 x 10"''/i^ Mpc'^ for 
the best-fitting HOD to the radio-detected LRG sample. 
This can now be compared to the measured density of 
2SLAQ LRGs which are currently detected as radio source^, 
0.041 ±0.003 X 10"* /i^ Mpc"^. Taking the ratio of these two 
numbers suggests that about 30% of all 2SLAQ LRGs with 
the same clustering properties as the radio-detected sam- 
ple are currently radio-loud. This in turn implies a fairly 
high duty cycle for radio-loud AGN in the central galax- 
ies of clusters at z ~0.55, and is consistent with the radio 
detection rate of the brightest 2SLAQ LRGs (28±8% for 
Mo.2,r < -23.75 from Table 4 of Sadler et al. 2007). 

The above discussion relies on our assumption that the 
form of the HOD for the radio-detected LRGs is the same as 
the standard form for a luminosity-limited sample. In that 
case the only mass-dependent probability of an LRG being 
radio-loud is determined by Mmin , and there is then a halo 
mass-independent probability of ~30% of a galaxy being 
currently powerful enough to be detected as a radio source at 
> 10^'' WHz"\ However it may be the case that the HOD 
for radio LRGs takes a different form from that given in 
equations l9land [T0l since one might expect it to consist of the 
standard form plus some term describing how the fraction 
of radio loud LRGs depends on halo mass (-F'rad(M)). 

We can attempt to determine _F'rad(M) by trying to 
modify the best fitting HOD to the matched sample in such 
away that we reproduce the clustering and space density 
of the radio-detected LRG population. We chose to model 
^rad(A^) as a power law such that 

F,ad(M) =/,(Af/M,ad)'', (19) 

where M^ad = 10"/i"^Mq and Frad(A/) ^ 1. 

Figure [7] shows the best fitting power law for -F'rad(Af), 
HOD and space density distribution, along with the best fit- 
ting standard HOD given in Table]!] The fit to the clustering 
is shown as the dotted line in Figure [S] The best fit has a 



^ This value is found simply by dividing the total number of 
radio-detected LRGs by the survey completeness and co-moving 
volume between 0.44 < z < 0.76 in the 2SLAQ survey area. Our 
clustering analysis uses the full 2SLAQ LRG sample, which is 
magnitude— limited and so includes a smaller contribution from 
lower-luminosity galaxies than would be expected in a volume- 
limited sample. For this reason, we argue that the simple space 
density derived above is the appropriate one to use in this compar- 
ison. If we integrate over the radio luminosity function (RLE) in 
Table 5 of Sadler et al. (2007), which uses the Vmax estimator to 
correct for the effects of the magnitude limit, we obtain a slightly 
higher value of 0.067 ± 0.004 X IQ-'^h^ Mpc'^ for the volume- 
limited sample, which would increase the estimated radio-loud 
fraction from ~ 30% to ~ 45% 
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M,„„(h-'M^) 

Figure 7. The mean number of LRGs per halo (top), the mean 
number of LRGs per halo normalised by the number of haloes 
(middle) , the fraction of radio— loud LRGs (bottom) , all as a func- 
tion of halo mass. The solid line shows the best fit to the matched 
sample, the dashed line shows the best fit to the radio sample 
when the density is matched, and the dotted line shows the best 
fitting matched HOD plus power-law radio fraction to the radio 
clustering. The red area shows the Icr error on the dotted line. 

slope (3 = 0.65±0.23 and normalisation fr = 0.14±0.02. The 
best fitting bias and effective halo mass are almost identical 
to those determined with the simple luminosity limit HOD, 
with bg = 2.93 ± 0.19 and Me// = 10.3 ± 1.9 x10^^/i"^Mq. 
This is reassuring, since it implies that the determination of 
the typical halo mass is largely independent of the form of 
the HOD, providing that it yields a reasonable fit to the 
clustering and space density. 

We also show in the bottom panel of Figure [7] -Frad(M) 
determined by dividing the best-fitting luminosity limit 
HODs to the radio and matched samples (shown by the 
blue and red lines in Figure [6}. While both HOD fits, with 
and without the power law, provide good statistical fits to 
our data, they diverge for halo masses above 10^^ H'^Mq 
(where the power-law fit suggests that ~60% of central 
cluster galaxies would be radio-loud, while the luminosity- 
limit HOD suggests that the fraction levels off at ~ 30%) 
and below a few times W^"^ Mq (where the luminosity- 
limit HOD suggests that there should be a sharp cut off in 
the fraction of radio-loud LRGs). Since the current 2SLAQ 
LRG data sample contains very few objects with halo masses 
above 10^^ H'^Mq or below 3 x 10^^/i"^Mq, a larger sample 
of radio galaxies is needed to test between different models 
in this regime. 



Table 5. Radio-detected fraction of brightest— cluster galaxies 
(above a limiting 1.4 GHz radio power Pi. 4) in the local universe 
and at z ^ 0.55. 



Sample 


Pi. 4 (WHz-i) 


Reference 




> 1023 


> 1024 




BCGs, z < O.ll 


33.3% 


22.2% 


Burns 1990 


BCGs, z < 0.2 


32.7% 


19.9% 


Lin & Mohr 2007 


BCGs, SDSS 


~ 30% 


- 20% 


Best et al. 2007 


2SLAQ LRGs 




> 30% 


This paper 



6 DISCUSSION AND CONCLUSIONS 

We have measured the two-point cross-correlation function 
of a sample of radio-detected LRGs at 2: ~ 0.55 and a 
sample of radio-quiet LRGs that matches the luminosity, 
colour and redshift distribution of the radio loud sample. 
We find that radio-detected LRGs at 2 ~ 0.55 are signifi- 
cantly more clustered than the matched radio-quiet sample, 
with clustering scale lengths (ro) of 12.3±1.2 and 9.02±0.52 

Mpc respectively. This result suggests that the radio- 
detected LRGs typically occupy more massive haloes than 
other LRGs of the same optical luminosity and stellar mass. 
We confirm this by fitting HODs to the observed cluster- 
ing, and show that the radio-detected LRGs have a typical 
halo mass of 10.1±1.4 xIO^/i'^Mq and a bias of 2.96±0.17, 
compared to a halo mass of 6.44±0.32 X 10"/i^^MQ and a 
bias of 2.49±0.02 for the radio-quiet LRGs. 

The clustering of the radio-detected LRGs is best fitted 
by a HOD containing only galaxies which occupy the cen- 
tres of haloes and are not satellites, though the current data 
do not allow us to exclude models with some radio satellite 
galaxies. We show that the dependence on radio fraction on 
halo mass can be modelled as a power-law with slope 0.65 
± 0.23, although it will take a larger samples of radio galax- 
ies with a more precise clustering measurement to rule out 
other forms of this dependence. Such samples are available 
at lower redshift, such as the SDSS/NVSS/FIRST sample 
presented by Best et al. (2005), which would also allow a 
measurement of any evolution in this relationship. Mandel- 
baum et al. (2008) use both a clustering and weak lensing 
analysis of this sample to show that the the halo masses of 
radio loud galaxies are about twice the mass of a compara- 
ble radio quite sample, in good agreement with the results 
presented here. 

Table [5] compares the radio detection rates for the 
most optically-luminous brightest cluster galaxies (BCGs) 
in the local universe with the 2SLAQ LRG value derived 
in §5.4. The values found in the three local studies are re- 
markably consistent, and imply that the radio power above 
which at least 30% of BCGs are detected rises from about 
10^^ WHz-i at z ~ 0.1 to at least 10^* WHz'^ at z ~ 0.55. 

As noted by Johnston et al. (2008) low-power (< 
10^^ WHz~^) radio galaxies in the 2SLAQ LRG sample have 
stellar populations which are generally no different from 
those in a radio-quiet comparison sample, and the major- 
ity of these radio galaxies have weak or no emission lines in 
their optical spectra. These properties are consistent with 
a population of 'low-excitation' radio galaxies powered by 
Bondi accretion of hot. X-ray emitting gas from the inter- 
galactic medium (Allen et al. 2006; Hardcastle, Evans & 
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Croston 2007). If this is correct, then the highly-clustered 
environment of the 2SLAQ radio galaxies is not surprising. 
The radio emission from central cluster galaxies could be be 
enhanced by the presence of a denser intracluster medium 
(ICM) which would both confine the radio lobes and pro- 
vide a more efficient working surface for the radio jets, thus 
boosting the observed radio luminosity (Barthel & Arnaud 
1996; Allen et al. 2006). However, Best et al. (2006) argue 
that this boosting of radio emission from cluster galaxies 
only applies to powerful radio sources which extend well be- 
yond the host galaxy and into the intracluster medium. Since 
most of the radio-detected 2SLAQ LRGs are low-power ra- 
dio galaxies whose radio emission is generally compact and 
confined within the host galaxy, this ICM-related boosting 
may not play a strong role. 

The most likely alternative is that the radio emission 
from central cluster galaxies is boosted because of a higher 
fuelling rate of gas onto the central black hole. This is plau- 
sible because gas from both the galaxy's own X-ray halo 
and a larger-scale cluster cooling flow can contribute to the 
fuelling of an active nucleus in central cluster galaxies (Best 
et al. 2007). 
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APPENDIX A: FITTING THE HALO MODEL 
TO A MAGNITUDE LIMITED SAMPLE 

In section [5] we estimate the bias and effective halo mass of 
effectively magnitude-limited galaxy samples using the halo 
model framework. The halo occupation distribution we use 
is designed to apply to volume limited galaxy samples and so 
its application here may result in biased estimates of these 
parameters. In this appendix we use the millennium simula- 
tion to investigate the magnitude of this bias and show that 
it is likely to be only ~5%, and comparable to the likely 
systematic error introduced by any differences between the 
assumed form of the HOD and the actual form. 

We construct two samples of galaxies generated using 
the latest version of the Durham group semi-analytic galaxy 
formation model detailed in Font et al. (2008). For the first 
sample we simply select all galaxies in the z — 0.02 volume 
with SDSS r magnitudes < —22 i.e. a volume limited sample 
of galaxies with a space density equivalent to our full LRG 
sample. We also wish to generate a magnitude limited sam- 
ple, and so calculate the apparent SDSS r magnitude of each 
galaxy in our volume assuming an observer at one edge of 
that volume. We then limit this sample to r < 19.67 which 
matches the space density of the volume limited sample. 
Both samples contain ~ 20, 000 galaxies. Fig. I All shows the 
absolute magnitude distribution and real space 2pt-auto cor- 
relation function for these samples. As one would expect the 
magnitude limited sample which contains some intrinsically 
fainter galaxies shows a slightly lower clustering amplitude. 
In order to test whether we are able to accurately measure 
the effective halo mass by fitting a HOD to the clustering 
measurements we first determine the actual HOD of the two 
galaxy samples shown as the solid lines in Fig. IA2I We then 
take a suitable analytic form of the HOD and fit it to both 
the actual HOD from the simulation and the derived 2pt- 
correlation function and space density, shown as the dashed 
and dotted lines respectively. It is worth noting that the ac- 
tual HOD derived from GALFORM does not exactly match 
our analytic form at low masses and would not match any 
form typically used in the literature. The clear 'bump' at 



Table Al. The Effective halo mass in units of IQ^^h-^MQ. 



Method 


Volume Limited 


Magnitude Limited 


Actual 


6.63 


6.38 


HOD fit 


6.98 ± 0.09 


6.87 ± 0.08 


Clustering fit 


7.01 ± 0.11 


6.78 ± 0.11 



; Vol limit HOD .: 




Figure A2. The measured HOD for the volume (red) and mag- 
nitude (blue) limited samples, as well as analytic HOD fits to the 
measured HODs (dashed lines) and the clustering (dotted lines). 



10 < Mhalo < 10 is caused by the onset of AGN feed- 
back in model (Carlton Baugh, Private Communication). It 
will be interesting to see if there is any evidence for this in 
directly determined HODs of massive galaxies. 

Despite the differences in the analytic and actual HODs 
the analytic HOD does provide good fits to the clustering 
and space density of the samples within the errors. Table 
I All gives the values M^ff determined from fitting the ana- 
lytic HOD to the measured HOD and the clustering as well 
as that determined directly from the simulation. M^ff is 
overestimated by ~7% fitting the analytic form to either 
measured HOD or clustering for both the volume and mag- 
nitude limited samples, due to the excess of galaxies in the 
'bump' in the HOD at low halo masses. The difference be- 
tween the fits to the HOD and the clustering is larger for the 
magnitude-limited sample, but is still only at the 2% level 
and is consistent with the error on the fit to the clustering. 
We can therefore be confident that the likely magnitude of 
any systematic error on Me/ / caused by fitting the HOD to 
a magnitude limited rather than volume-limited sample is 
less than the measurement errors for all the samples con- 
sidered in Section [S] It is interesting to note that a larger 
systematic error in M^ff would be introduced if the true 
HOD is not well represented by the assumed analytic form 
of the HOD. However, one would expect still relative mea- 
surements of Meff to be representative even in this case. 



12 



Wake et al. 




Mr (Mag) r (h ^ Mpc) 

Figure Al. The absolute magnitude distribution (left) and the real space 2pt-corrclation function (right) for the volume limited (dashed 
line) and magnitude— limited (solid line) simulated galaxy samples. 
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